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First-principles calculations with the full-potential linearized augmented plane-wave (FP-LAPW) 
method have been performed to investigate detailed electronic and linear optical properties of 
Sr2Cu02Cl2, which is a classical low-dimensional antiferromagnet (AFM) charge transfer (CT ) 
insulator. Within the local-spin-density approximation (LSDA) plus the on-site Coulomb interac- 
tion U (LADA+(7) added on Cu Sd orbitals, our calculated band gap and spin moments are well 
consistent with the experimental and other theoretical values. The energy dispersion relation agrees 
well with the angle resolved photoemission measurements. Its linear optical properties are calculated 
within the electric-dipole approximation. The absorption spectrum is found to agree well with the 
experimental result. 

PACS numbers: 71.15.Ap, 74.25.Jb, 74.25. Gz 
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I. INTRODUCTION 

The study of the strongly correlated electronic sys- 
tems in two dimensions (2D) has been a very active and 
interesting problem since the discovery of the high 
superconductors, in which calculation of the electronic 
structures of these materials by the first principles meth- 
ods is a first step toward a microscopic understanding 
of them. Among these materials, the Sr2Cu02Cl2 is an 
idea paradigm since it has an extremely weak coupling 
between its Cu02 planes, and is thus a more exaggerated 
quasi-2D system than the most other high-Tc materials. 

Experimental and theoretical efforts have been made to 
investigate electronic structures of the Sr2Cu02Cl2. For 
example. Wells et al.^ made a measurement on the angle- 
resolved photoemission spectroscopy (ARPES) above its 
Neel temperature 256.5 ± 1.5K and got its energy band 
dispersion relations, which are in consistent with the 
t — J model calculations along some k-paths but not 
all. To achieve more precise results, Kim et alwk man- 
aged successfully to do the ARPES experiments on the 
Sr2Cu02Cl2 at 150K, which is below its Neel tempera- 
ture, and found that it would be better to include the 
second and third nearest neighbor hopping terms, i.e., to 
use the t — t' — t" — J mode&iiSiSii, for well describing the 
observed phenomenon. The first principles calculation is 
an important theoretical method to get information on 
the ground state of materials. Hua Wu et al& used the 
linear combination of atomic orbital (LCAO) method to 
calculate the electronic structure of the Sr2Cu02Cl2 in 
an AFM structure. Their results give out reasonable den- 
sity of states (DOS) to some extent, but do not show the 
dispersion relationship. D. L. Novikov et al^ also stud- 
ied the Sr2Cu02Cl2 by FP-LMTO numerical calculation 
without taking into account its magnetic structure and 
the strong correlation effect of Cu 3d electrons. 

So, in this paper we try to do a more detailed study 
on the AFM Sr2Cu02Cl2 with the FP-LAPW^" method 



within the local-spin-density approximation and plus the 
on-site Coulomb interaction on Cu 3d orbitals. Then, 
based on this and within the electric-dipole approxima- 
tion, its linear optical responses are also calculated. A lot 
of worksiiii^iiii^ii^, mostly are model calculations and 
experimental measurements, had been done to study its 
optical properties for many years in order to understand 
the microscopic mechanism of the high-Tc superconduc- 
tivity exhibited in the layered copper oxides. As to our 
best knowledge, there is still no first-principles calcula- 
tion study on its optical properties. Our result gives out 
a broad peak at the charge transfer (CT) gap, which is 
the dominant optical absorption feature in these mate- 
rials, and other peaks in higher energy range, which are 
also well comparable with the experimental data. 



II. METHODOLOGY 

FP-LAPW~ method treats all electrons without the 
shape approximations for the potential and charge den- 
sity. The Perdew- Wang's exchange-correlation energies 
are used in the local density approximation (LDA)iS,. In 
our calculation we have taken the crystal and magnetic 
structures of the Sr2Cu02Cl2 from Ref. [17]. As shown 
in Fig. 1, in the antiferromagnctic phase, the unit cell 
is doubled while the first Brillouin zone (BZ) shrinks to 
half. The radii of the atomic spheres are taken as 2.0 a.u. 
for Cu, CI, and Sr and 1.6 a.u. for O, respectively. And 



the cutoff of the plane- wave basis is set to RmtKn 



=7.0 



(the Kraax is the plane-wave cut-off, and the Rmt is the 
smallest radius for all atomic sphere radii.). Thus, the 
resulting number of the plane waves used in our calcula- 
tion is about 4950 plus 190 local orbitals. Self-consistent 
calculations are performed with 200 k-points in the first 
BZ by the tetrahedron method. For checking the numeri- 
cal convergence, we increase k-points to 500 and find that 
the difference between the two converged total energies 
is about 0.02 mRy. In order to include the strong corre- 
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FIG. 1: The crystal and magnetic structure of the AFM 
Sr2Cu02Cl2 material. Cul means the Cu atom with spin 
up and Cu2 spin down. Also shown is the two-dimensional 
Brillouin zone. Solid (dashed) line is the one without (with) 
consideration of the AFM magnetic structure. 



lation effect of tlie Cu 3d electrons, we use tlie LSDA+C/ 
approachl^ with U=6.Q eV and J=125 meV. And the 
case of U=7.5 eV is also calculated for comparison. Here 
the U,J values are taken from Ref. [8], [19] and [20]. 
Convergence is assumed when the difference between the 
input and output charge densities is less than 0.0001 
e/(a.u)'^ and the difference of total energies in the last 
two iterations is less than 0.0001 Ry. Optical responses 
are calculated within the electric-dipole approximation. 
The imaginary part of the dielectric function can be ex- 
pressed by 

e2H-:^EEI<c,fc|e.p|t;,fc>i2 



c,v k 

E,{k)-niu], 



where c and v represent the conduction and valence 
bands, respectively, and \c,k >, \v,k > are the eigen- 
states obtained from the FP-LAPW calculations, p is 
the momentum operator, e is the external electric field 
vector, and lu is the frequency of incident photons. We 
still use 200 sampling k-points in the irreducible wedge 
and the linear tetrahedron scheme improved by Blochl 
et al^ for the BZ integrations. Using Kramers-Kronig 
(K-K) transformation, we can get the real part of the 
dielectric function, and further calculate the absorption 
spectrum of the Sr2Cu02Cl2. 



III. RESULTS AND DISCUSSION 

The calculated electronic structures for the AFM 
Sr2Cu02Cl2 with t/=6.0 cV are plotted in Fig. 2, 3 and 
4. From Fig. 2, it is seen that there is an indirect band 
gap of about 1.2 eV with the top of the valence band at 
M(7r/2, 7r/2) point and the bottom of conduction band 
at X(7r, 0) point, which is smaller than the experimental 
value of 1.9 eV^, but much closer to it than the value of 
0.84 eV when U= 5 eV in Ref. [8]. If a larger f/=7.5 
eV is selected and other parameters are unchanged, an 
indirect gap of 1.68 eV is obtained, which is very close 
to the value of 1.63 eV obtained by the LCAO method-. 
The energy dispersion relation agrees very well with the 
results in Ref. [1] and [2]. But, along the direction from 
X(7r, 0) to r(0, 0) it is better than that obtained from 
the t — J model calculationi. The energy dispersion of 
the top-most valence band along r(0, 0) to M(7r/2, 7r/2) 
reaches to about 0.45 eV, larger than the experimental 
value of 280±60 meVi But calculation with C/=7.5 eV 
shows that a larger U will suppress the dispersion to 0.32 
eV, which is much closer to the experimental value. Anti- 
ferromagnetic order of Sr2Cu02Cl2 can be clearly seen 
from Fig. 3(a) and 3(b). The total spin is zero because 
both of Cul and Cu2 have anti-parallel spin moments of 
about 0.59 per atom, which is also bigger than the 
experimental value 0.34 ^g^. It is found that the larger 
U of 7.5 eV will increase the spin moment of Cu atom to 
0.63 per atom. 

Although the hybridization of O 2p and Cu ^dx2_y2 
between -3 and -1.2 eV can be very clearly seen in Fig. 
3(b) and 3(c), it is still obvious that the topmost valence 
band is mostly composed of the O 2p orbital, and a care- 
ful analysis shows that near the top of the valence band 
the O 2p component is about 62.8%. Oppositely it is 
seen from Fig. 4(a) that the Cu 3^3.2 _j^2 occupies more 
of the lowest conduction band with its component being 
about 69.4%. In the case of U=7.5 eV, the O 2p com- 
ponent increases to 65.6% at the top of the valence band 
and Cu 3d^2_y2 occupies more to 70.6% of the lowest 
conduction band, which indicates that Sr2Cu02Cl2 is a 
charge transfer insulator. Thus the transition from O 2p 
to Cu 3dx2-y2 orbitals will dominate the absorption edge. 
The CI 3p orbital hybridizes with Cu 3d orbital in the 
range of from -5 to -3 eV, mostly with Cu dyz + d^z and 
d^z^-r^ orbitals (see Fig. 3(d), Fig. 4 (b) and (d)), all 
which are in the plane perpendicular to the Cu-0 plane. 
The hybridization between CI and Cu is less than that of 
Cu and O because O 2p orbitals extends more than CI 
3p, and so hybridizes more with Cu 3d orbitals. 

Study of optical property can give more informations 
on the electronic structures. In Fig. 5, absorption spec- 
trum calculated with ?7=6.0 eV are presented. Although 
U—7.5 eV gives the similar spectrum shapes, but all the 
peaks would shift to higher energy by about 0.4 eV. So, 
we prefer to take U=6.0 eV as an example in the follow- 
ing discussion. Comparing it with the experimental data 
given by R. Lovenich et.alM- and Moskvin et.alM', we 
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FIG. 3: Total and projected density of states (DOS) of AFM 
Sr2Cu02Cl2. Vertical dashed line refers to the Fermi level, 
(a) Total DOS for spin up and down; (b) partial DOS of Cul 
3d; (c) and (d) are DOS of O 2p and CI 3p, respectively. 
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FIG. 2: The band structure of the insulating AFM 
Sr2Cu02Cl2. The Fermi energy is set to zero. Z is (0, 0, 
tt) and other high symmetry k-points are defined in Fig. 1. 



find that they match very well near the charge transfer 
gap of about 2eV, corresponding to electrons transferred 
from O 2p to Cu 3d orbitals, which can be clearly seen 
from the band structure in Fig. 2 along the path from 
X(7r, 0) to M(7r/2, 7r/2). Around 2.5 eV, there is one 
lower peak in both our calculation and Ref. [15], which 
was deemed by Moskvin et. alMi as the other part of the 
double-peak structure of the lowest excitation. The third 
strong and broad peak appears at about 3.5 eV. Accord- 
ing to the idea of Pothuizcn et alm^, a nonbonding band 
(NBB), mostly O 2p^ orbitals, lies at about 1.5 eV be- 
low the Zhang-Rice band (ZRB), so the energy needed 
to transfer an electron from NBB to Cu 3d is about 3.5 
eV, which is also shown by Moskvin et al. in their ex- 
perimental observations. Again, the broad peak around 
4.6 eV appears both in our fist-principles result and their 
measurements. Since the Cu 'idj.2_y2 occupies the lowest 
part of the conduction band from Fermi level to about 
1.5 eV, and CI 3p dominates the most states between -4.5 
and -2.5 eV, the quite broad peak around 4.6 eV can be 
ascribed to the hopping of electrons from CI 3p orbitals 
to Cu 3dj.2_y2 orbitals, which can be clearly seen from 
Fig. 3 (b), (d) and Fig. 4 (a). Usually, the CUO4 model 



Fig. 4 



2 
1 

-1 

-2 
2 
1 

> 

0) 

12 -1 

^ -2 
S 3 

^ 2 

O 1 
Q 
-1 
-2 
-8 
4 
2 

-2 
-4 
-6 





(a) d , , 


" 1 i 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 " 


1 °l' -<)!'«i||j|iyl^^^ 1 1 1 






1 1 1 1 1 1 1 
1 1 1 1 1 1 1 " 


■ 1 1 1 l| ;l ' 1 1 1 1 1 1 1 1 1 1 


(c)d 


' 1 1 1 Ml 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 






1 1 1 1 1 ' 1 ' 


1 1 1 1 1 1 1 



-7 -6 -5-4-3-2-101234 

Energy (eV) 



FIG. 4: Spin polarized projected DOS of Cul and Cu2 3d 
electrons. Solid line represents that of Cul and dotted line, 
Cu2. Also the vertical dashed line refers to Fermi level, (a) 

d^2_y2, (b) d^^2_^2, (c) d^y, {d)dyz + d^^^. 
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FIG. 5: The calculated absorption spectra of AFM 
SraCuOaCla. 

calculation on the Sr2Cu02Cl2 has just neglected the CI 
effect by replacing the CI atom by the O atom, while 
here we have shown the contribution of CI atoms to the 
optical properties of the Sr2Cu02Cl2. 

Finally, we want to emphasize that correct choice of U 
value is critical for the LSDA+?/ calculation. One crite- 
rion is whether the chosen U value can give a reasonable 
physical result in consistent with the experimental obser- 



vation. Combining our obtained results with that in Ref. 
[8] for U—5.0 eV, we can see that U=6.0 and 7.5 eV are 
both physically reasonable although the former has been 
selected for the main discussion on Sr2Cu02Cl2's optical 
properties. And it is clear that larger U influences Cu 3d 
orbital more than O 2p, and causes a wider gap, less band 
dispersion and a bigger magnetic moment on Cu. Also, 
larger U shifts the absorption edge to higher energy. 



IV. CONCLUSION 

In summary, we have given a more detailed first- 
principle calculations on the band structure and linear 
optical properties of the AFM Sr2Cu02Cl2 by using the 
LSDA-f [/ approach in the FP-LAPW formahsm. The 
obtained results indicate that it is a charge-transfer insu- 
lator with an indirect gap of about 1.2 eV, which is less 
than the experimental value. The E vs. k relationship 
agrees very well to the ARPES observations. Also, we 
have found that our calculated linear optical susceptibil- 
ity is well consistent with the experimental data not only 
around the charge-transfer gap, but also in the higher 
energy range. 
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